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Abstract. We consider the problem of verifying the existence of ground states of the 
ID nonhnear Schrodinger equation for an interface of two periodic structures: 

^u" + V{x)u- Xu^T{x)\u\P-^u onR 

with V{x) = Viix),T{x) = ri(a;) for a; > and V{x) = V2{x),T{x) = T2{x) for x < 0. Here 
Vi,y2,ri,r2 are periodic, A < mincr(-^ + V), and p > I. The article [T. Dohnal, M. 
Plum and W. Reichel, "Surface Gap Soliton Ground States for the Nonlinear Schrodinger 
Equation." Comm. Math. Phys. 308, 511-542 (2011)] provides in the ID case an existence 
criterion in the form of an integral inequality involving the linear potentials Vi , V2 and the 
Bloch waves of the operators — ^-I-Vl,2 — A. We choose here the classes of piecewise constant 
and piecewise linear potentials Vi^2 and check this criterion for a set of parameter values. In 
the piecewise constant case the Bloch waves are calculated explicitly and in the piecewise 
linear case verified enclosures of the Bloch waves are computed numerically. The integrals 
in the criterion are evaluated via interval arithmetic so that rigorous existence statements 
are produced. Examples of interfaces supporting ground states are reported including such, 
for which ground state existence follows for all periodic ri_2 with esssupri_2 > 0. 



1. Introduction 

An interface between two nonlinear periodic media in the n— dimensional nonlinear Schrodinger 
model can act as a waveguide so that localized solutions, so called surface gap solitons (SGS), 
exist as shown analytically in [3]. Experimentally such waveguiding has been demonstrated 
in nonlinear photonic crystals, see e.g. [9|[TT|[T2]. There are also a number of numerical obser- 
vations of SGSs in the ID and 2D nonlinear Schrodinger equation (NLS), see e.g. [HEllSlE] . 

In [3] an existence criterion for strong ground states of the n— dimensional NLS 

(n-NLS) (-A + V{x) - \)u = T{x)\u\p~\, x 

was proved with V{x) = Vi{x),T{x) = Ti{x) for Xi > and V{x) = V2{x),T{x) = T2{x) for 
xi < under the condition A < m.ma{—A + V). The functions Vi,V2,Ti,T2 are assumed 
periodic in each coordinate direction and the exponent p satisfies p E (1, 2*), where 2* = 
if n > 3 and 2* = cxDif?T, = l,2. A strong ground state is defined to be a minimizer of 
the corresponding total energy J^„ |(|Vup + V{x)u'^) — -^^lul^^^dx restricted to the Nehari 
manifold N = {u E H\W) \ {0} : \Vu\^ + {V{x) - X)u^ - T{x)\u\P+^dx = 0}. The 



Date: February 17, 2012. 

2000 Mathematics Subject Classification. Primary: 35Q55, 78M30, 65G20; Secondary: 35J20, 35Q60. 
Key words and phrases, nonlinear Schrodinger equation, surface gap soliton, ground state, variational 
methods, interface, periodic material, verified numerical enclosures. 

1 



2 



TOMAS DOHNAL, KAORI NAGATOU, MICHAEL PLUM, AND WOLFGANG REICHEL 



results of [3] include sufficient conditions for the existence of strong ground states. These 
conditions involve information about the strong ground stated wi,W2 of the purely periodic 
problems (n-NLS) with y = Vi, F = Fi on and = V2,F = F2 on respectively. 
In the case n = 1 these conditions could be formulated in terms of the Bloch waves of the 
two purely periodic linear problems. As neither the ground states wi,W2 nor the Bloch 
waves are generally known explicitly, [3] did not produce explicit examples of ground state 
supporting interfaces except for an example where the potentials are related by scaling: 
Vi{x) = k'^V2{kx), Ti{x) = 7^F2(fcx) with certain conditions on k and 7, see Theorem 5 
in [3]. All other existence examples were asymptotic; either in A or in Fi — F2. 

The most practical existence criteria in [3] are those for the ID case n = 1. In this article we 
provide a number of explicit ID examples of interfaces satisfying these criteria. We consider, 
therefore 

(1.1) - u" + V{x)u -\u = T{x)\u\P-\ on R 

with 



(1.2) V{x) 
and 

(1.3) F(x) 



Vi(x), X > 0, 

V2{x), X <0, 

Ti{x), X > 0, 

F2(x), x<0 



under the assumptions 

(HI) Vi,V2,Fi,F2 are 1-periodic, 
(H2) ess sup Fi > 0, 2 = 1,2, 
(H3) 1 < p < 00, 
(H4) X<mina{-£ + V), 

which were needed in [3]. 

Next, recall the criterion given in Theorem 7 in [3] for the existence of SGS ground states 

of dUD. 

Theorem 1. Assume (H1)-(H4) and for i = 1,2 define by q the energy of a strong ground 
state of (- + Vi{x) - X)u = Ti{x) \u\p~'^u on R. 

(a) // ci < C2, then a sufficient condition for the existence of a strong ground state of 
ffLTj) IS 



;i.4) h := j (V2{x) - Vi{x)^u^]:\xf dx < 0, 



where u_ (x) = p_ {x)e'^^^ , with ki > and pi^ periodic, is the Bloch mode decaying 



at -00 of + Vi{x) - A. 



^The existence of strong ground states of the purely periodic problem on E" was proved in [5]. 
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(b) // ci > C2, then a sufficient condition for the existence of a strong ground state of 



where (x) = {x)e '^'^^ , with K2 > andp\_ periodic, is the Bloch mode decaying 
at +00 of-^ + V2{x)-X. 

When the ordering of ci, C2 is unknown, Theorem [T] can still be used by establishing nega- 
tivity of both of the integrals Ji and l2'- 

Corollary 2. // both /i,/2 < 0, then a strong ground state exists irrespectively of the order 
of ci and 02, and thus, of the choice o/ri,r2 and p (within the assumptions (H1)-(H3)). 

As seen from (11. 4p and (II. 5p . any ordering Vi{x) > V2{x) or Vi{x) < V2{x) on [0, 1] leads 
to Ji < < /2 or /2 < < Ji, whence the assumption of corollary [2] is not satisfied. 

If information on the ordering of the ground state energies ci,C2 is available, the corre- 
sponding criterion (a) or (b) in Theorem [1] can be checked. This is the case, for instance, 
with the dislocation interface 



where Vq, Fq are 1— periodic, ess sup Fq > 0, and where Ti,T2 G [0,d] are the dislocation 
parameters. In this case Ci = C2 so that we have 

Corollary 3. For the interface (II. 6p suppose Ii < or I2 < 0. Then there exists a strong 
ground state irrespectively of the choice of p (within (H3)) and ofV^ (within the above as- 
sumptions). 

We will use direct constructive approaches to verify the respective conditions (II. 4p and 
(II. 5p . These require mainly the Bloch waves of the two purely periodic linear problems on 
M. We consider two types of potentials Vi^2'- piecewise constant and piecewise linear. For 
piecewise constant potentials we calculate the needed Bloch modes in closed form by hand, 
so that we can check these conditions directly. For piecewise linear potentials we use a 
computer-assisted approach, i.e. we compute verified enclosures for the Bloch modes, and 
use these to enclose the integrals Ii and I2 in (ll.4p and (II. 5p . These computer-assisted results 
are completely verified and thus give a rigorous mathematical proof since all numerical errors 
are taken into account. In principle, even much more general potentials can be treated by 
this approach; we have chosen piecewise linear ones for simplicity. 

In the rest of the paper we impose the condition 



ffLTjl IS 



(1.5) 




(1.6) 




A < min{inf Vi, inf V2} 



which ensures (H4) without having to actually calculate the spectrum. 

In Section [2] below we consider interfaces with piecewise constant potentials Vi and V2 and 
in Section [3] we study the case of piecewise hnear potentials. 
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2. PiECEWiSE Constant Potentials Vi and V2 

When the potentials Vi and V2 are piecewise constant, the integrals /i, I2 can be calculated 
explicitly although the closed form involves the inverse of a transcendental function. We 
calculate the formulas for /i, I2 explicitly and evaluate these numerically for a set of parameter 
values. The evaluation is done in interval arithmetic (using the Matlab toolbox Intlab [ID]). 
The resulting values of /i, I2 are thus enclosed in intervals and when the supremum of such 
an interval is negative, the corresponding integral Ji or I2 is then verified to be negative. 



2.1. Bloch Waves for a Piecewise Constant Potential. 

Let 

fa, < X < s, 
^°W = U, s~<.<l 

with a,b eR and s G (0, 1). 

The Bloch waves of —u" + Vo{x)u = Am on R have the form 

u±{x) = p±{x)e^'^^ with p±{x + 1) = p±{x) , k > 

since A lies in the resolvent set of — ^ + Vq, see At the same time, due to the piecewise 
constant nature of Vn 



(2.2) u±{x) 



^±gxv^ + ^^e"^^^, < X < s. 



^±gxv^ + ^±e-^^^, s < a; < 1. 



The vectors are determined via the condition for u±{x) at x = s and the condition 
that the Floquet multipliers of u± are e^'^ respectively. For we thus obtain the system 



/a~\ e»^°^ -y^a~\ e-"^^^ -y/b^ ^sVb^^ v^F^ g-a^b^ 



— e 



Solving det(y4(K)) = yields 

cosh(K) = — (Va — A + Vfe — A)^ cosh (s\/ a — A + (1 — s)-\/fe — A 

Ay a — Xyh — A L V 

— (V a — A — — A)^ cosh ^sV a — A — (1 — s)^/h — A 



(2.3) 



INTERFACES FOR SGS EXISTENCE IN ID NLS 



The solution vector ^+ is proportional to 

p—s\/ a— A 

^1+ 



(2.4) 



sinh((l - s)^/b^) 



= {a-b) 



QS\/b—X—K Q\/b—X~s\/a~X ' 



A + Va - Wb - A 



^(V a— A— A)s— K g— A 

p(\/a^^+\/5— A)s— re p^/b—X ' 



^4*" = A — a + — AVfe — A. 
The system for reads y4(— k),^" = 0, so that 
(2.5) ^-(k)=^+(_;,). 

2.2. Dislocation Interface. 

Let us consider the dislocation interface fll.6p with ti = T2 =: t for the piecewise constant 
potential 

f a, < X < 1/2, 
6, 1/2 < X < 1 

with a, 6 G M. In this case u^^\x) = u±{x + r) and ■uli^^(x) = u±{x — r) with m± in (12.21) . ^'^ 
in ([23D, ([23]) and k in (O, where we set s = 1/2. 
A direct integration then produces for < r < 1/2 

b — a 



-2k 



2Va^ 



^-^g(l+2r)Va-A ( ^{2r-l)^ _ 1 ) ( 1 _ 



-2T%/a-X 



(2.6) 



+ 



b — a 



-2k 



2Vb^ 



^-^^il+2r)Vb^fil-2r)Vb^_^] ( ^2rVb^X _ ^ 



+ ^fe-^'+^^^^ ( e(2-i)v^ - 1 ) ( 1 



and for 1/2 < r < 1 



b — a 



-2k 



(2.7) 



+ 



2Va^ 



b — a 



^-2g(l+2r)v/^ /'g(l-2T)V^ _ /'g2(T-l)V^ _ I 



i-A 



-2k 



2Vb^ 



^3 e- 



2(r+l)VF:A /'g(l-2r)v^ _ ;l ] ( 1 _ g 



2{t-1)v/F^ 
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For the integral I2 we have for < r < 1 /2 
a — h 



h{r) 



^+2g(l-2r)v'^ /g(2T-l)v^^ _ ^ 



(2. 



+2^(2r-l)V^ / „(1-2t)V^ 



2V6^ 



+2 {2r-l)v^ /^{2r-l)v^ 



and for 1/2 < r < 1 



/2M 



a — b 



(2.9) 



+ 



2Va - A 



a — b 



2Vb^ L 



^+2g(l-2r)v^ ^g(l-2r)V^ _ ;L 
_^ ^+2^(2r-l)v^^ ^^{2t-1)V^ _ 
"^+2g2(l-r)v^ j^g(l-2r)v^ _ I 



1 - e 



-2rVa-A 



g2rv'a-A _ ^ 



2rVb^ 



1 - e 



-2TVb~X 



2{T-l)Va^ _ ^ 
I _ g2(l-r)v^^ 
_ g2(r-l)v^ 
,2(l-r)v^ir:A _ 



For the dislocation interfaces with a = 1 and b G {2, 6} Fig. [T] shows regions of the (r, A) 
plane where the integral Ii or I2 is negative, i.e. where ground state existence is guaranteed. 
These regions were computed using interval arithmetic. The domain [0, 1] x [—2, 0.98] in 
the (r, A)-plane was completely coverecfl by two dimensional intervals (squares) of size 1 / 600 
along each dimension and when for a given square the interval arithmetic evaluation produced 
Ik < (which means a negative supremum of the enclosure of Ik), the square was shaded. 
Note that it is a priori clear that both Ii and I2 are zero at r = 0, 1/2, and 1 because for 
these values Vo(x + r) — Vo{x — r) = due to the 1— periodicity of Vq and the integrands in 
/i, I2 thus vanish . The use of interval arithmetic in the computations then necessarily results 
in small neighborhoods of r = 0, 1/2, 1 where the sign of the integrals cannot be determined. 

As Fig. [T] shows, ground state existence is guaranteed in both cases 6 = 2 and 6 = 6 in 
almost the entire parameter domain (r. A) G [0,1] x [—2,0.98]. Interestingly, in both cases 
we find regions where Ji and I2 are negative. In these regions ground state existence thus 
follows by Corollary [2] irrespectively of the nonlinear coefficients Fi and r2 and of the order p 
because both criteria in Theorem [T] are satisfied. The coefficients and the order can be freely 
changed (within (H1)-(H3)) and the ground state will still exist! In particular, ri,r2 do not 
have to be related by the dislocation property in (11. 6p . 

In Fig. [2] and [3] the enclosures of \/a — \\/b — A Ji and \l a — \\lb — A I2 as functions of r 
are plotted for the two examples in Fig. [T]at two values of A. The quantities /i,/2 contain 
ya — A and a/6 — A in the denominator. To reduce the amount of round off error for A close 
to min{a,6}, we plot \/ a — \\Jb — A I\^2- Note that in Fig. [2] (b) the r-region, where both 



^AU intervals whose endpoints are not floating-point numbers are safely enclosed in slightly larger intervals 
by use of Intlab. 
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(a) a=1,b=2 (b) a=1,b=6 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 1. Regions in the (r, A) plane, where suprema of the enclosures of Ii 
and I2 are negative for the dislocation interface with a piecewise constant Vq 
and a = 1. (a) 6 = 2; (b) 6 = 6 (computations performed in interval arith- 
metic). 

Ji and I2 are negative, is larger than at the corresponding value of A in Fig. [1] (a). This is 
because Fig. [T]is produced using two dimensional intervals (intervals in the (r, A)-plane) for 
input while in Fig. |2]the input are one dimensional intervals. The larger overestimation due 
to interval arithmetic in the case of two dimensional intervals leads to a smaller region where 
negativity of Ji and I2 is verified. 

In Fig. [3] (b) a small left neighborhood of r = 1/2 is visible where the infima of the 
enclosures of both Ji and I2 are positive so that existence of strong ground states cannot be 
concluded based on Theorem [H 



(a) a=l, b=2, X=0 ,g-4 (b) a=l, b=2, X=0.9i 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

T T 



Figure 2. Scaled plots of the interval enclosures of \/ a — \\/h — A Ji and 
yja — \\/b — A I2 as functions of r corresponding to Fig. [1] (a). In (a) A = 
and in (b) A = 0.94 (computations performed in interval arithmetic with the 
r-interval width 1/600). 
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-1 
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-4 
-5 
-6 
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(a) a=l, b=6, X=0 



(b) a=l, b=6, X=0.94 



0.2 



0.4 0.6 

T 



0.8 




Figure 3. Scaled plots of the interval enclosures of Va — Xy/b — X Ii and 
\/ a — X\/h — X I2 as functions of r with a and b as in Fig. [1] (b). In (a) A = 
and in (b) A = 0.94 (computations performed in interval arithmetic with the 
r-interval width 1/600). 



2.3. Interface with no Knowledge about Ci,C2. 



Here we consider a general interface (11. 2p . (II. 3p with the piecewise constant structure (12. ip 
for both Vi and V2. In detail, Vi is given by (12. ip with (a, 6, s) replaced by (ai,fei,Si) and 
V2 is given by (12. ip with {a,b,s) replaced by (02,^2, ■§2). For simphcity we choose the jump 
locations in the middle of the periodicity cell: Si = S2 = 1/2. 

The Bloch waves are now given by (12. 2p and (l2.3p -f E3]) with (a, 6, s) replaced by 

(ai, 61, 1/2). We denote the resulting k in (12. 3p by ki. Analogously we obtain u^f^ and denote 
the resulting vectors in (l2.4p .( E3]) by C^. 

Because the ordering of Ci and C2 is unknown in this case, Theorem [T] can be used to prove 
ground state existence only if both Ji and I2 are negative. If this occurs, the existence of 
a strong ground state is then completely independent of the nonlinear periodic coefficients 
Fi, r2 and of p (within (H1)-(H3)). We show below that similarly to the dislocation example 
such cases can be found here as well. 

The integrals from Theorem [T] now become 



(2.10) 

= e-2-i(a2-ai)ere2" + 



+ e-'^'{b2-b,)^^^^ + 



e ^'"^(aa - aQ 
2Vai - A 
(62 - 61) 



^-^(e^_l)+^-^(l-e- 



-2/ti 



2y%^ _ ^Vb^^ 



+ ^4 (e 
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and 



(2.11) 



h = (ai - a2)Ci+C2^ + 




where is the same as in (12. 4p with (a, b, k, s) replaced by (02, &2, ^2, 1/2), k,2 is the same 
as K in (12. 3 p with (a, 6, s) replaced by (02, 62, 1/2). 

Fig. H] shows regions of the (ai, A) plane where the integrals Ji, I2 are negative. The shaded 
region is where both Ji and I2 are negative, i.e. where ground state existence is guaranteed 
irrespectively of the coefficients ri,r2 and of p (within (H1)-(H3)). Similarly to Fig. [T]we 
covered the region {(ai,A) G : — 1 < ai < 1.5,-1.5 < A < min{ai, 02, ^1, ^2} — 0.0256} 
completely with squares of size 0.0056 in each dimension and used interval arithmetic to 
evaluate Ji and I2. 



Figure 4. Region in the (ai,A) plane, where suprema of the enclosures of 
/i, I2 are negative for the interface (II. 2p with piecewise constant Vi and V2 and 
bi = 2.5, 02 = 2, 62 = 2, si = S2 = I (computations performed in interval 
arithmetic). 



In Fig. [5] the enclosures of a/c^i ~ AV^i — A Ji and y/a2 — \ Vb2 — A I2 as functions of ai 
are plotted for the example in Fig. |l]at two values of A. At A = 0.8 the value of I2 is always 
positive while at A = both Ji and I2 are verified negative for ai G (0.11, 0.85), where ground 
state existence thus follows. 



b, =2.5, 82=1, b2=2 




-0.5 







0.5 



1.5 



3. Piecewise Linear Potentials Vi and V2 



Here we consider continuous piecewise linear functions as potentials. Unlike in the case of 
piecewise constant potentials in Section [21 explicit formulas for the Bloch modes u± are now 
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0.04 r 

0.02 

- 

-0.02 
-0.04 
-0.06 



(a) 


b = 
1 


2 


.5, a^^ 




J a — 


T 


/b^Ii 




^\ 


/b 




ai = 


= A 







X 10" 



(b) b =2.5, a=l, b =2, X=0 . 




1 0.01 Vqi - AVti - a /i 

v'q2 - AV'?2 - A /2 

■ = A 




Figure 5. Scaled plots of the interval enclosures of i/ai — X\/bi — X Ii and 
\/a2 — X\/b2 — A /2 as functions of ai with 02,^1 and 62 as in Fig. HJ In (a) 
A = and in (b) A = 0.8 (computations performed in interval arithmetic with 
the ai-interval width 1/1000). 

generally not available. We compute the Bloch waves via the numerical enclosure method 
presented in [7]. All presented results are therefore verified. 



Example 1. Let us consider the dislocation interface (11. 6p with ti = T2 =: t for the following 
1-periodic potential 

f 4x- 1, xE [0.25,0.5], 
Vo{x) = l -4x + 3, x G [0.5,0.75], 

[0, xe [0,0.25] U [0.75,1], 

cf. Fig. El 




Figure 6. Potential Vo{x) on its period for example [T] 



Let M+ and U- be the fundamental solutions of 

-u" + {Vo{x) -X)u = 
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with the form 

(3.1) u+{x) = e~'^''p+{x), u-{x) = e'^''p-{x), 

where k > is the characteristic exponent and p+,P- are 1— periodic functions. Note that 
for p := > 1 

u+{x + 1) = p~'^u+{x), U-{x + 1) = pu-{x) 

hold. 

Since Vo{x) is an even function, u+{—x) is also a fundamental solution. x) must be 

a linear combination first of both U-{x) and but since grows at — oo whereas 

x) and U-{x) both decay at — oo, the factor of in the linear combination must be 

zero, i.e. x) = cu-{x) holds for some c G C. Noting that the fundamental solutions can 
be normalized, we can define 

:= x) 

after computing U-{x). Then by simple calculations we see that Ji = I2 holds. 

In Fig. [7] the enclosure of Ji as a function of r is plotted for the two cases: (a) A = — 1, 
(b) A = —0.1. In case of (a) Ji is negative for r G [0.04, 0.49] and in case of (b) Ji is negative 
for r G [0.2,0.49]. Hence ground states exist in these cases. 




0.2 0.4 0.6 0.8 1 ■ 0.2 0.4 0.6 0.8 1 



Figure 7. Plot of the interval enclosure of Ji as a function of r with interval 
width 0.01 for example [T] 

Example 2. As another example of the dislocation interface (11. 6p with ti = T2 =: t we 
consider the 1-periodic potential 

( 30x -4.5, X G [0.15,0.25], 

Vo{x) = I -30x + 10.5, X G [0.25,0.35], 

[0, X G [0,0.15] U [0.35,1] 

shown in Fig. [HI 
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0.15 0.35 0.5 0.65 0.85 1 



Figure 8. Potential Vo(x) on its period for example [2] 



Because Vq is not even, Ji is generally different from I2 and we plot in Fig. [9] the enclosure 
of both Ji, I2 as functions of r for the case (a) A = — 1 and (b) A = —0.01. Several regions are 
observed, where Ji or I2 are negative and where by Theorem [T] ground state existence follows 
for the corresponding interfaces fll.6p . As in Section [2] we also find intervals where both Ii 
and I2 are negative so that ground states exist for arbitrary Fi, F2 and p (within H1-H3). In 
case (a) both Ji and I2 are negative for r G [0.01,0.31] U [0.51,0.66] and in case of (b) both 
Ji and I2 are negative for r G [0.18,0.31]. 




Example 3. Here we consider a general interface fll.2|] with piecewise linear Vi,V2. As 
explained in Section 12. 3[ we need to show that both Ii and I2 are negative in order for 
Theorem [T] to yield ground state existence. As we noted below Corollary [21 we should 
violate a monotone order between Vi and V2 in order to possibly obtain negative Ii and I2 
simultaneously. 
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(3.2) V,{x) 

(3.3) V2{x) 
see Fig. [IO](a). 



2000X-400, a; G [0.2,0.25], 
-2000a; + 600, x e [0.25,0.3], 
0, a; G [0,0.2] U [0.3,1], 

20X-14, X G [0.7,0.75], 
-20a; + 16, x G [0.75,0.8], 
0, X G [0,0.7] U [0.8,1], 



100 
80 
60 
401- 

20 








(a) and 



■IOV2 



0.2 



).4 0.6 0.8 

X 



-0 .01 



-0.03 



(b) and 



Figure 10. (a) Potentials Vi{x) and V2{x) in fj3l2|) . fl3l3|) on their period, 
(b) Interval enclosure of /i,/2 as functions of A with interval width 0.01 for 
potentials in (O, (Q. 



We have verified that for A G [—2, —0.01] both Ii and I2 are negative. We plot in Fig. [10] 
(b) the enclosure of both Ji, I2 as functions of A using intervals with width 0.01. 

In case that supp(V2) is closer to supp(Vi) as follows (see Fig. [TT]), we verified that Ji 
is negative and I2 is positive for A G [—2, —0.01], whence ground state existence cannot be 
concluded using Theorem [1] 

( 2000X-400, X G [0.2,0.25], 

(3.4) Vi{x) = I -2000X + 600, x G [0.25, 0.3], 

[0, X G [0,0.2] U [0.3,1]. 

( 20X-6.5, X G [0.325,0.375], 

(3.5) V2{x) = I -20X + 8.5, x G [0.375, 0.425], 

[0, X G [0,0.325] U [0.425,1]. 
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Figure 11. Potentials Vi{x) and V2{x) in (13. 4 p and (13. 5p on their period. 
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